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Summary. Uniaxially anisotropic antiferromagnets in a field along the easy axis are 
studied with the help of ground state considerations and Monte Carlo simulations. 
For classical models, the XXZ model as well as variants, we analyze the role of non- 
collinear spin configurations of biconical or bidirectional type interpolating between 
the well-known antiferromagnetic and spin-flop structures. Possible experimental 
applications to layered cuprate antiferromagnets are discussed. Finally, results of 
quantum Monte Carlo simulations for the 5=1/2 XXZ model on a square lattice are 
presented, and compared with previous findings. 

1 Introduction 

Since many years uniaxially anisotropic antiferromagnets in a field have been 
studied extensively, both experimentally and theoretically. The magnets are 
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Fig. 1. Spin orientations on neighboring sites showing AF (a), SF (b,d), and BC 
(c) as well as BD (e) ground state structures in XXZ (a,b,c) and anisotropic XY 
(a,d,e) antiferromagnets. 
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known to display at low temperatures, T, upon increasing the field H along 
the easy axis, antiferromagnetic (AF) and spin-flop (SF) phases [1]. A proto- 
typical model describing these phases is the XXZ Heisenberg antiferromagnct, 
with the Hamiltonian 

wxxz = j J2 [ A ( s ? s ? + s i sv j) + s i s i ] - H Y, s i w 

where Sf, Sf , and Sf denote the spin components at lattice site i. The first 
sum runs over pairs of neighboring sites (i, j) of the square or cubic lattice; 
J > is the exchange integral, and A, < A < 1, determines the strength 
of the anisotropy along the easy axis (z-axis) . The field H acts along the z- 
axis. Classical XXZ antiferromagncts on square and cubic lattices have been 
analyzed using Monte Carlo techniques first about three decades ago [2, 3]. 

Additional phases in the (H, T)-plane, observed in experiments and theo- 
retical studies, have been attributed to, for instance, further anisotropy terms 
and/or interactions ranging beyond nearest neighbors [4, 5, 6]. 

One of the main aims of the present contribution is to draw attention 
to recent theoretical analyses of the classical XXZ antiferromagnet and its 
analogue for spins with only two components, the anisotropic XY antiferro- 
magnet [7, 8, 9, 10, 11]. Especially, the importance of non-collinear struc- 
tures of biconical (BC) or bidirectional (BD) type, see Fig. 1, is emphasized 
[10, 11, 12]. By adding a single-ion anisotropy term to the XXZ model, 
these structures may be enhanced or suppressed, depending on whether 
that term favors a planar or a uniaxial anisotropy [11]. Some of the recent 
theoretical analyses have been partly motivated by experiments on quasi- 
twodimcnsional cuprate antiferromagncts, the 'telephone number compounds' 
(Ca, LajuCuziO^i [13, 14, 15, 16, 17]. Thence, we shall discuss also more 
complicated models for uniaxially anisotropic twodimensional antiferromag- 
nets proposed to describe such compounds, in particular Ca 5 LagCu2iO 41 
(here, one may mention previous and recent experimental studies on related 
quasi-twodimensional antiferromagncts [18, 19, 20, 21, 22] as well). 

Finally, we shall consider the quantum, 5* = 1/2, version of the XXZ 
antiferromagnet, which is equivalent to a Bose-Hubbard model [23, 24], on 
a square lattice, being, especially, of current interest in the context of super- 
solids. The model is simulated using the method of stochastic series expansions 
(SSE). New results [10] on the phase diagram will be compared with previous 
ones [24, 25, 26]. 

Our findings will be summarized at the end of this contribution. 

2 The classical XXZ Heisenberg antiferromagnet 

The classical XXZ Heisenberg antiferromagnet is described by the Hamilto- 
nian given in Eq. (1), where we shall deal with spin vectors of length unity on 
square and cubic lattices. 
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Fig. 2. Phase diagram of the classical XXZ model with A= 2/3 (a) and 4/5 (b). 



We first consider the two-dimensional version. Phase diagrams in the 
(T, -ff)-plane are depicted in Fig. 2, where we set the exchange anisotropy 
A to be equal to 2/3 and 4/5, the latter case being the standard choice 
[2, 3, 7, 8, 10]. 

The general topology of the phase diagram seems to be independent of the 
concrete value of A, < A < 1, comprising the long-range ordered AF and the 
algebraically ordered SF phase. The boundary lines to the disordered phase are 
in the Ising universality class for the AF case, and in the Kosterlitz-Thouless 
universality class [27] for the SF case [2, 7]. The AF and SF phase boundary 
lines approach each other very closely near the maximum of the SF phase 
boundary in the (T, 7?)-plane, see Fig. 3. Accordingly, at low temperatures, 
there may be either a direct transition between the AF and SF phases, or two 
separate transitions with an extremely narrow intervening phase may occur. 

Indeed, recent simulations suggest a narrow (disordered) phase between 
the AF and SF phases [7] , extending presumably down to zero temperature [8] . 
The evidence has been provided by determining the universality classes of the 
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Fig. 3. Phase diagram of the classical XXZ model near the maximum of the bound- 
ary line of the AF phase, A = 4/5. 

transitions at low temperatures [7, 8] , remaining to be of Ising- or Kosterlitz- 
Thouless type, and by finite size arguments in the limit of T approaching zero 
temperature [8]. The presence of that intervening phase may be argued [10] 
to be closely related to the highly degenerate ground state occurring at the 
field H c i = 4J\/I — A 2 separating the AF and SF structures. At that point, 
not only AF and SF configurations have the same energy, but also biconical 
structures. Those structures may be described by the tilt angles Oa and Ob 
characterizing the orientations of spin vectors at neighboring sites, belonging 
to the two different sublattices A and B, see Fig. 1. The two tilt angles are 
interrelated by [10, 11] 



with the BC configurations interpolating continuously between the AF and 
SF structures, where the tilt angle Oa ranges from to ir. 

The relevance of BC fluctuations in the transition region between the AF 
and SF phases at low temperatures may be conveniently seen by studying 
probability functions of the tilt angles, such as the probability p2 (Oa i Ob) for 
finding the two angles, 0a and Ob, at neighboring sites and the probability 
p(0) for encountering the tilt angle O [10, 11, 12]. An illustration is depicted 
in Fig. 4, showing that the line of local maxima in p2 follows closely Eq. (2), 
signaling that the degenerate BC structures are present in that region. The 
rather low probability for configurations deviating only slightly from the AF 
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structure, see Fig. 4, is caused by the small radii of their 'cones'. In fact, at 
T = 0, the probability of P2 along the line given by Eq.(3) is proportional to 
sj sin{0 A)sin(0 b) ■ In any event, the, presumably, narrow disordered phase 
seems to be governed by degenerate BC fluctuations, with a hidden 'tetracrit- 
ical point' at zero temperature [10]. 
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Fig. 4. Probability p2{0A, Ob) showing the correlations between the tilt angles Oa 
and Ob on neighboring sites for the XXZ antiferromagnet, with 80 x 80 spins, 
at H/J = 2.41, fcsT/J = 0.255, and A = |. p2 is proportional to the gray scale 
[10]. The superimposed black line depicts the relation between the two angles Oa 
and Ob in the ground state, Eq. (2). 

For the classical XXZ antiferromagnet on a cubic lattice early renormal- 
ization group arguments [28, 29] and Monte Carlo simulations [3] suggested 
that the triple point, at which the AF, SF and paramagnetic phases meet, is 
a bicritical point with 0(3) symmetry (obviously, such a point is excluded to 
occur, at T > 0, in two dimensions due to the well-known theorem by Mermin 
and Wagner [30]). Only a few years ago, this scenario has been questioned, 
based on high-order perturbative renormalization group calculations [31]. It 
has been predicted that, instead of the bicritical point, there may be a 'tetra- 
critical biconical' [29] point, due to an intervening ordered 'biconicaP phase 
in between the AF and SF phases, or a point at which first-order transition 
lines meet. 

Our previous Monte Carlo simulations [7] for the thrcedimcnsional XXZ 
antiferromagnet agreed with a first-order transition between the AF and SF 
phases at low temperatures. Based on our current simulations [32], again 
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for A = 0.8, we locate a triple point at k B T t /J = 1.025 ± 0.015 and 
H t / J = 3.88 ± 0.03, in reasonable agreement with the estimate by Landau 
and Binder [3]. So far we did not specify its (multicritical) character. More- 
over, we analyzed p2 1 showing that BC fluctuations occur in the transition 
region between the AF and SF phases temperatures well below T t . However, 
in contrast to the twodimensional case, one now observes a tendency towards 
coexistence of the AF and SF phases, as reflected by fairly pronounced max- 
ima at the corresponding points in the {0a, Ob) plane, compare with Fig. 
4. We also estimated critical exponents from monitoring the size dependence 
of the maxima in the longitudinal, i.e. along the easy axis, as well as the 
transverse staggered susceptibilities and the specific heat near the transition 
between the AF and SF phases at temperatures somewhat below that triple 
point. In agreement with previous findings [3, 7] and the behavior of p 2l a 
transition of first order between the AF and SF phases is strongly suggested 
[32]. 

3 Variants and applications to layered cuprate magnets 
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Fig. 5. Phase diagram of the anisotropic XY antiferromagnet with A — 0.8. 

In the following, classical variants of the twodimensional XXZ antiferro- 
magnet are studied. From a theoretical point of view, one may like to check 
the robustness of the topology of the phase diagram against modifying the 
model, to identify genuine features. From an experimental point of view, one 
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may like to have a (semi-)quantitative description of measurements. Address- 
ing the first aspect, we added a single-ion anisotropy term to the XXZ model, 
and we also analyzed the anisotropic XY antiferromagnet on a square lattice. 
Addressing the experimental aspect, we focused on layered cuprate antiferro- 
magnets, the socallcd 'telephone number compounds' (Ca, La) 14CW24O41, in 
particular Ca^La^Cu240n. 

3.1 Classical anisotropic XY antiferromagnet 

Reducing the number of spin components to two and keeping the uniaxial 
anisotropy, one arrives at the anisotropic XY antiferromagnet, with the Hamil- 
tonian 

^xy - j [ s ? s * + AS f s i ] - H E 5 * ( 3 ) 

where the x-axis is now the easy axis. As before, < A < 1. We set A = 0.8. 
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Fig. 6. Probability P2(Oa,Ob) for the anisotropic XY antiferromagnet with A= 
0.8 for a system with 100 x 100 lattice sites in the transition region between the AF 
and SF phases at ksT/J — 0.558 and H/J = 2.44. P2(Oa,0b) is proportional to 
the grayscale. The superimposed solid line depicts the relation between the two tilt 
angles Oa and Ob in the ground state, see Eq. (2). 

The topology of the phase diagram looks like in the XXZ case [11, 12], 
compare Figs. 2 and 3 with Fig. 5. The AF and SF phase boundary lines 
approach each other very closely near the maximum of the AF phase boundary 
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in the (T, iJ)-plane. Accordingly, at low temperatures, there seems to be, 
again, either a direct transition between the AF and SF phases, or two separate 
transitions with an extremely narrow intervening phase. 

In principle, now a bicritical point, with 0(2) symmetry, may occur at 
non-zero temperature, being of Kosterlitz-Thouless type. Our simulations, 
however, suggest that, like in the XXZ case, there is a narrow disordered 
phase intervening between the AF and SF phases down to temperatures well 
below the point where the AF and SF phases approach each other very closely 
[11, 12]. In particular, critical exponents of the staggered susceptibilities are 
found to be compatible with the Ising universality class holding now for the 
transitions of the AF as well as the SF phases to the disordered phase. The 
narrow intervening phase seems to be due to degenerate bidirectional struc- 
tures, see Fig. 1, arising from the highly degenerate ground state at the field 
H c i, which separates the AF and SF structures. This behavior is in complete 
analogy to the one for the XXZ antiferromagnet. 

At the highly degenerate ground state, the two tilt angles ©a and 0b 
of the BD structures are interrelated analogously to Eq. (2). In contrast to 
the XXZ antiferromagnet, however, the probability P2(Oa, Ob) along the line 
of the interrelated tilt angles is now constant at T = 0. Both properties, 
the degeneracy and the (almost) constant value of P2, tend to hold at low 
temperatures as well. This is displayed in Fig. 6, depicting p 2 in the transition 
region between the AF and SF phases at low temperatures. 

3.2 Adding a single ion anisotropy to the XXZ model 

The classical XXZ model, Eq. (1), on a square lattice is modified by adding a 
single-ion anisotropy term of the form 



which either, D < 0, enhances the uniaxial anisotropy A, or, D > 0, may 
weaken it due to a competing planar anisotropy. The sign of D will have 
drastic consequences for the phase diagram [11, 12, 33]. Due to the single- 
ion term, the highly degenerate ground state, at H c i in the XXZ model, is 
removed, suppressing altogether the BC structures, when D < 0, or spreading 
them over a finite range of fields, limited by H c \ a and H c ib, when D > 0, see 



In the latter case of a competing anisotropy, an ordered BC phase arises 
at low temperatures, bordered by the AF and SF phases, as shown in Fig. 
7. Based on renormalization group calculations [34, 35, 36, 37], the transition 
between the BC and SF phases may be argued to be in the Ising universality 
class, while the transition between the BC and AF phases is expected to be 
in the XY universality class, being the Kosterlitz-Thouless universality class 
in two dimensions. This description is in accordance with our simulational 
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Fig. 7. Phase diagram of the XXZ antiferromagnet with a competing single-ion 
anisotropy, A = 0.8 and D/J — 0.2. 




Fig. 8. Histograms for the probability of the tilt angle p(0) for the XXZ antiferro- 
magnet with a competing single-ion anisotropy, D/J = 0.2, at ksT / J — 0.2, at the 
fields given in the inset. Lattices with 80 x 80 spins are simulated. 
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data, as inferred from critical exponents for staggered susceptibilities and 
magnetizations at the two different phase boundary lines [11, 12, 33]. 

In the BC phase the interrelated tilt angles are changing continuously, 
at fixed low temperature, with the field. This behavior is displayed by the 
probability function p(0), as illustrated in Fig. 8. By increasing the field, 
the peak positions correspond first to the AF structure, shifting gradually 
towards each other, reflecting BC structures, and finally merging in one peak 
characterizing the SF phase. 

As seen in Fig. 7, the extent of the BC phase shrinks with increasing 
temperature. Eventually, the BC phase may terminate at a tetracritical point 
[5, 34, 35, 36, 37], where the AF, SF, BC, and paramagnetic phases meet [11]. 




k D T/J 



Fig. 9. Phase diagram of the XXZ antiferromagnet with a single-ion anisotropy 
fostering the uniaxiality, A = 0.8 and D/J = —0.2. 

In the case of a negative single-ion anisotropy, D < 0, enhancing the 
exchange anisotropy, there are no ground states of BC type. In Fig. 9, a 
typical phase diagram is depicted, where D/J = —0.2, showing long-range 
ordered AF and algebraically ordered SF phases. 

At low temperatures, we observe a transition of first order separating the 
AF and SF phases. Evidence for that kind of transition is provided, especially, 
by the critical exponent describing the size-dependence of the maximum in 
the longitudinal staggered susceptibility and by a coexistence phenomenon of 
AF and SF structures in the transition region between the two ordered phases, 
showing up, e.g., in p(Q) [11, 12]. 



Uniaxially anisotropic antiferromagnets 11 



When the uniaxiality is solely due to a single-ion anisotropy, A = 1, D < 0, 
BC structures do not occur as ground state. Accordingly, one may tend to 
believe that, at low temperatures, a direct transition of first order between 
the AF and SF phases takes place, in contrast to a conflicting claim [38]. Of 
course, this aspect needs to be clarified. 

Here, attention is also drawn to interesting recent work on a twodimen- 
sional Heisenberg antiferromagnet with long-range dipolar interactions pro- 
viding a uniaxial anisotropy [39]. 



3.3 Descriptions related to Ca^LaQCu^O^x 



The quasi-twodimensional uniaxially anisotropic Heisenberg antiferromagnet 
Ca^La<jCui^O^\ shows intriguing magnetic features, as the consequence of 
an interplay of spin and charge properties in the coupled CuO^ spin chains 
[13, 40], see below. Perhaps most interestingly, at low temperatures a sharp 
transition at a fairly low field along the easy axis is followed, at a higher field, 
by an anomaly in the susceptibility. 




To describe the measured spin-wave dispersion of Ca$Laa,Cu2±0^\ Mat- 
suda et al. proposed a classical twodimensional model with short-range ex- 
change interactions and a single-ion anisotropy [15]. However, the model does 
not reproduce thermal properties of that magnet like the sharp transition 
[16, 42]. That transition has been argued to indicate the thermal breaking of 
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'defect stripes' [13, 40, 41]. The 'defects' correspond to nonmagnetic Cu 3+ - 
ions, due to mobile holes, replacing some of the magnetic Cu 2+ -ions in the 
Cu0 2 spin chains (the defect concentration is about 10 percent). Indeed, such 
a transition has been described by an Ising model with mobile defects [14]. 
Introducing nonmagnetic, mobile defects in the model of Matsuda et al. is, 
however, not sufficient to reconcile the discrepancy with the measurements 



The experimentally observed anomaly at higher fields has been explained 
qualitatively as indicating the onset of merely local spin-flop structures re- 
lated to a significant decrease in the mobility of the defects or holes [40]. 
This presumption has been used in a twodimcnsional uniaxially anisotropic 
Heisenberg antiferromagnet, with an exchange anisotropy and short-range 
competing interactions by including defects quenched at randomly chosen lat- 
tice sites [17]. The model parameters have been carefully chosen, partly on 
theoretical, partly on experimental grounds. In fact, the model then repro- 
duces (semi-)quantitatively the field dependence of the anomaly when chang- 
ing the temperature, see Fig. 9. From the simulations, one easily sees that the 
anomaly is driven by the onset of merely local spin-flop structures, its local 
character being due to the quenched random defects. 

Note that the competing exchange interactions may induce helical spin 
configurations when tuning the parameters suitably [17]. 

We should like to mention a recent study on these 'telephone number 
compounds' applying density functional theory [44] , which might also be useful 
to quantify model parameters. 

4 S = 1/2 XXZ quantum antiferromagnet on a square 



Previous Monte Carlo simulations of the S=l/2 XXZ antiferromagnet on a 
square lattice suggest that there is, at low temperatures, a direct transition 
of first order between the AF and SF phases [24, 25, 26]. 

Schmid et al. [24] performed quantum Monte Carlo simulations to de- 
termine the phase diagram, for A = 2/3. They found a topology which re- 
sembles that of the classical XXZ antiferromagnet with a negative single-ion 
anisotropy, compare with Fig. 9. On the boundary line of the AF phase a 
tricritical point has been reported to occur at fcsT tc / J w 0.141, below which 
the transition to the paramagnetic phase is of first order. The triple point, at 
which the AF, SF, and disordered phases meet, is proposed to be a critical 
endpoint, located at fcsT ce / J w 0.118. 

To check these predictions, we performed large scale quantum Monte Carlo 
simulations [10] using the method of stochastic series expansions with directed 
loop updates [45]. From those simulations, considering larger system sizes 
and improved statistics, we conclude that the previous analysis has to be 
viewed with care. For instance, we studied the model, A = 2/3, at fc^T/ J= 



[43]. 
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Fig. 11. Positions of the maxima of the magnetization histograms as a function 
of the inverse system size, 1/L, for the S=l/2 XXZ antiferromagnet on a square 
lattice, with A = 2/3. The inset exemplifies two histograms for systems of linear 
size L = 32 (circles) and L = 150 (squares) at fesT/J = 0.13 and the coexistence 
fields H/J = 1.23075 and H/J = 1.232245 [10]. 



0.13, i.e. in between T tc and T ce , near the AF phase boundary. In particular, 
we monitored the size dependence of peak positions in the magnetization 
histograms, see Fig. 11. Obviously, the two peaks, corresponding to AF and 
SF structures, may well coincide at the transition in the thermodynamic limit. 
Thence the transition may well be continuous, in contrast to the previous 
suggestion. Actually, at the lowest temperature we studied, fc^T/J w 0.096, 
a continuous transition might still occur [10]. We conclude that the previous 
[24] scenario with the triple point being a critical endpoint needs to be shifted 
to somewhat lower temperatures, if it exits at all. 

In any event, a clue on possibly distinct phase diagrams for the classical 
and quantum XXZ antiferromagnets on a square lattice may be the possibly 
different role of biconical fluctuations. That aspect deserves further studies. 



5 Summary 

In this contribution we presented results of recent Monte Carlo simulations 
on classical XXZ antiferromagnets in a field along the easy axis as well as 
classical variants and on the S=l/2 XXZ antiferromagnet on a square lattice. 

Basic aspects of phase diagrams and applications to Ca^La^Cu^Ow are 
discussed. 
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The role of non-collincar structures of biconical and bidirectional type 
in classical models is emphasized. These structures have an important effect 
on phase diagrams, in particular, the transition region between the AF and 
SF phases at low temperatures, and they may provide a clue to explain the 
possibly different topology of the phase diagrams of classical and quantum, 
S=l/2, XXZ antiferromagnets in two dimensions. 

We thank especially A. Aharony, K. Binder, B. Biichner, B. Kastcning, R. 
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for cooperation, help, remarks, information, and/or encouragement. Financial 
suppport by the Deutsche Forschungsgemeinschaft under grant SE324/4 and 
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